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Summary. We consider asymmetric (nonreciprocal) wave transmission through a 
layered nonlinear, non mirror-symmetric system described by the one-dimensional 
Discrete Nonlinear Schrodinger equation with spatially varying coefficients embed- 
ded in an otherwise linear lattice. Focusing on the simplest case of two nonlinear 
sites (the dimer), we compute exact scattering solutions such that waves with the 
same frequency and incident amplitude impinging from left and right directions have 
different transmission coefficients. The stability of some particular solutions is ad- 
dressed. We show that oscillatory instability may lead to the formation of stable 
extended states coexisting with a nonlinear defect mode oscillating at a different 
frequency. Numerical simulations of wave packet scattering are presented. Gaus- 
sian wave packets with the same amplitude arriving from opposite directions on 
the dimer are indeed trasmitted very differently. Moreover, asymmetric transmis- 
sion is sensitively dependent on the input parameters, akin to the case of chaotic 
scattering. 



1 Introduction 

The possibility to control energy and/or mass flows using nonlinear features of phys- 
ical systems is a fascinating issue both from the point of view of basic science as well 
as from the applied one. In the context of wave propagation trough nonlinear media, 
the simplest form of control would be to devise a "wave diode" in which electro- 
magnetic or elastic waves are transmitted differently along two opposite propagation 
directions. 

In a linear, time-reversal symmetric system this possibility is forbidden by the 
reciprocity theorem. As stated by Lord Rayleigh in his treatise The theory of sound 

[1]: 

Let A and B be two points [...] between which are situated obstacles 
of any kind. Than a sound originating at A is perceived at B with the 
same intensity as that with which an equal sound originating at B would be 
perceived at A. In acoustics [...] in consequence of the not insignificant value 
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of the wavelength in comparison with the dimension of ordinary obstacles 
the reciprocal relation is of considerable interest. 

To achieve the desired effect one must thus violate the hypotesis of the theorem. 
In linear systems, a standard way is to break the time-reversal symmetry by applying 
a magnetic field as done, for instance, in the case of optical isolators. An entirely 
alternative possibility is instead to consider nonlinear media. At least in principle, 
this option would offer a whole new range of possibilities of propagation control 
based on intrinsic material properties rather than an external field. 

Asymmetric wave propagation induced by nonlinearity arises in several differ- 
ent domains. Among the first examples discussed in the literature is the asymmetric 
phonon transmission through a nonlinear interface layer between two very dissimilar 
crystals [2]. In the field of nonlinear optics a relatively vast number of approaches 
exist. A so-called all-optical diode has been proposed first in Ref. [3, 4] and later on 
realized experimentally [5]. There are also proposals to employ left-handed metama- 
terials [6], quasiperiodic systems [7], coupled nonlinear cavities [8] or VT— symmetric 
waveguides [9, 10]. Extension to the quantum regime in which few-photon states dis- 
play a diode effect has been proposed [11]. 

In the realm of acoustics the possibility of realizing a diode has been demon- 
strated for nonlinear phononic media [12, 13]. Another promising context is the prop- 
agation of acoustic pulses through granular systems. Indeed, experimental studies 
demonstrated a change of solitary wave reflectivity from the interface of two gran- 
ular media [14]. More recently, demonstration of rectification of mechanical energy 
at sonic frequencies in a one-dimensional array of particles has been also reported 
[15]. 

Despite the variety of physical contexts, the basic underlying rectification mech- 
anisms rely on nonlinear phenomena as, for instance, second-harmonic generation 
in photonic [16] or phononic crystals [12], or bifurcations [15]. In those examples the 
rectification depends on whether some harmonic (or subharmonic) of the fundamen- 
tal wave is transmitted or not. 

A related question is the possibility that the transmitted power at the same fre- 
quency and incident amplitude would be sensibly different in the two opposite prop- 
agation directions. In this Chapter we address the above problem with the Discrete 
Nonlinear Schrodinger (DNLS) equation [17, 18] with site-dependent coefficients. 
It has been demonstrated [19] that DNLS equation can be a sensible approxima- 
tion for the evolution of longitudinal Bloch waves in layered photonic or phononic 
crystals. Variable coefficients describe different nonlinear properties of each layer 
and the presence of defects. In the realm of the physics of cold atomic gases, the 
equation is an approximate semiclassical description of bosons trapped in periodic 
optical lattices (see e.g. Ref. [20] and references therein for a recent survey). Beyond 
its relevance in many different physical contexts, the DNLS equation has the big 
advantage of being among the simplest dynamical systems amenable to a complete 
theoretical analysis. For our purpose, it is particularly convenient as it allows to 
solve the scattering problem exactly without the complications of having to deal 
with wave harmonics [21]. 

In Sec. 2 we outline the model and show some examples of asymmetric plane- 
wave solutions. The issue of their stability is briefly addressed in Sec. 3. It is also 
shown that oscillatory instability may lead to the formation of stable extended states 
coexisting with a nonlinear defect mode oscillating at a different frequency. In Sec. 4 
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we report some numerical simulation of wave packet's scattering and illustrate its 
dependence on initial intensity. Finally, a brief summary of the results is given in 
Sec. 5. 



2 The model 

The DNLS equation with spatially varying coefficients, defined on an infinite one- 
dimensional lattice is given by 

i<j>„ = Vn<J>n ~ 4>n+l ~ 4>n-l + Ct„\(f> n \ 2 (f>„ , (1) 

where units have been chosen such as the coupling C = 1. 

We will assume the usual scattering setup where V„ and a n are non vanishing 
only for 1 < n < N. The two semi infinite portions (n < 1, n > N) of the lattice, 
model two leads where the wave can propagate freely [22] . Let us look for solutions 
of the associated stationary transmission problem 4> n (t) = u n exp(— ifxt) 



of the form 



/j,u„ = V„u„ — u n+ i — u n -i + a n \u n \ 2 u n 1 < n < N (2) 



R e lkn + Re- %hn n < 1 
Te ikn n>N 



where fi — — 2cosfe and < k < n for the wave coming from the left direction; 
Ro,R and T are the incident, reflected and transmitted amplitudes respectively. 
The solution sought must be complex in order to carry a non vanishing current 
J = 2|T| 2 sinfc. 

To break the mirror symmetry with respect to the center of the nonlinear 
portion, one must choose at least one of the two sets of coefficients V„, a„ such 
that V n 7^ Vn-u+i, ct-n otN-n+i- Note that the transmission of the right- 
incoming wave with the same _Ro and fx is computed by solving the problem with 
(Vn,a n ) — > (Vjv-n+i, ajv-n+i) (i.e. "flipping the sample"). In the following, we 
will adopt the convention to label with — k the right-incoming solutions with wave 
number k. Nonlinearity is essential as for a n = the transmission coefficient is 
the same for waves coming from the left or right side, independently on V„ due to 
time-reversal invariance of the underlying equations of motion [23]. 

The standard way to solve the problem is to introduce the (backward) transfer 
map [24, 25, 26, 22] 

Un-l = —V n + (V n - fl + a n \u n \ 2 )u n , V„-l = tin • (4) 

Note that these are complex quantities therefore the map is nominally four- 
dimensional. However, due to conservation of energy and norm, it can be reduced 
to a two-dimensional area-preserving map [24, 25, 26, 22] with an additional control 
parameter (the conserved current J). The solutions are straightforwardly found by 
iterating (4) from the initial point un = T exp(ikN) , vn = T exp(ik(N + 1)) dic- 
tated by the boundary conditions of Eq. (3). For fixed T and k, the incident and 
reflected amplitudes are determined as 

exp(-ifc)uo — vo exp(ik)u - v 

Hq — -. TTT , , . , K — 



exp(— ik) — exp(ifc) ' exp(ifc) — exp(— ik) 
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and the transmission coefficient is t{k,\T\ 2 ) = \T\ 2 /\Ro\ 2 Note that if («o,«o) = 
(un,vn) (periodic point of the map) then t = 1. 

For very short chains (oligomers), t can be computed analytically. For instance 
for the dimer N = 2: 



t = 



l + (6 2 



*)(e ik - <5i) 



(5) 



where 

For the trimer N = 3 
t = 

where 



<5 2 



V 2 -fi + a 2 T 2 , 

Vi - n + aiT 2 [l - 2<5cosfc + <$?]. 



e lfc - <5i + (e ifc - (5 3 )(1 - 62(61 - e tk )) 



(6) 



;fci2 



63 = V 3 -M + a 3 |T| 2 
<5 2 = V 2 -fi + a 2 \T\ 2 \S a -e 
61 = Vi - LL + ai \T\ 2 \l + 8 2 (e lk -S 3 )\ 2 . 

The formulas apply for k > (left-incoming waves); the transmission for right- 
incoming waves is obtained by exchanging the subscripts 1 and 2. 

Two examples of the dependence of t on the input power are shown in Fig. 1. 
The curves display a multistable behavior and, for strong enough intensities, are 
sizeably different indicating nonreciprocity. The effect is maximal in the vicinity 
of the nonlinear resonances that are detuned differently for the k > and k < 
cases yielding intervals of input values were multiple solutions exist only for one 
propagation direction [21]. The lower panel of Fig. 1 shows that even a moderate 
increase of the number of sites (N = 3) dramatically increase the complexity of the 
curves as expected due to the mixed phase-space of the underlying transfer map 
[24]. 

To conclude this section, we mention that an alternative approach would be 
to use the Green's function formalism previously used to compute the stationary 
states for an electron moving on a chain with nonlinear impurities [27]. Indeed, for 
the case of a symmetric nonlinear dimer, resonance phenomena are demonstrated 
that lead to complete transmission through the dimer [27]. Of course, we expect 
such an approach to yield the same results when applied to the present case. 



3 Stability of scattering solutions 

An important issue is the dynamical stability of the solutions. To the best of our 
knowledge no systematic study of scattering solutions of the type described above 
has been presented in the literature (see [28] for the case of NLS equation with 
concentrated nonlinearities continuum case and [29] for an analysis of a related 
problem, the nonlinear Fano effect). 

The linear stability analysis is performed [17, 18] by letting <p n — (u n + 
Xn) exp(— ifit), and linearizing the equation of motion to obtain 




Note that \ n is complex. Letting 

\„ = ^4n exp(iAt) + B* exp(-i\*t) 
then A is solution of the eigenvalue problem 

XA n = — (£n, — n)An + A n + 1 + An-l — 6nB n 

XB n — +(e n — A*)Bn — Bn+1 — Bn-l + #n^n (8) 

with £ n = V n + 2a n \u n \ 2 , 8 n = a n u^. Note that, at variance with the case of e.g. 
breather solutions, the solutions are complex, and also the coefficients 8 n in (8) 
are complex as well. As it is known, the eigenvalues come in quadruplets of the 
form ±A, ±A*. If eigenvalues have a nonzero imaginary part then the solution is 
unstable. Generally speaking, equilibria of an Hamiltonian system can lose spectral 
(and therefore linear) stability in two ways: a pair of real eigenvalues can either 



6 



Stefano Lepri and Giulio Casati 



(i) merge at the origin and split onto the imaginary axis (saddle-node bifurcation) 
or (ii) collide at a nonzero point and split off into the complex plane, forming a 
complex quadruplet (Krein bifurcation). The latter case correspond to an oscillatory 
instability. 

To solve the linear problem (8) exactly one should impose to the solution a def- 
inite plane wave form (with complex wave numbers) in the two seminiifinite linear 
parts of the chain. The matching of such waves through the nonlinear portion re- 
duces the infinite-dimensional problem (8) to an homogeneous linear system of 2N 
equations, whose solvability condition, along with the dispersion relations, yields a 
set of nonlinear equations for the unknowns. The details of this method (which is 
technically more involved than a straightforward diagonalization) will be presented 
elsewhere [30]. Here, we limit to illustrate the stability properties of some representa- 
tive solutions by a more direct approach, i.e. by solving numerically the eigenvalue 
problem for a finite truncation — N p < n < N + N p of the chain (M sites with 
M = 2N P + N ± 1), checking that the relevant eigenvalues of the resulting 2M x 2M 
matrix are not affected by the truncation error. 




-2-1012 -2-1012 
Re \ fte\ 



Fig. 2. Upper panels: real (squares) and imaginary (diamonds) parts of the two 
solutions marked with I and II in Fig. 1: they correspond to incident waves with 
the same input |-Ro| 2 = 2 having transmission coefficients t — 0.99 and t = 0.30 
respectively. Lower panels: the spectrum in the complex plane N p = 200. Isolated 
eigenvalues for for solution I are ±1.316±0.502i for II, ±2.60 and ±2.46 respectively. 



Of course, for small enough nonlinearities/amplitudes the solutions should be 
stable. Since our main object of interest here is in the large asymmetry effects, we 
concentrate on the cases of strongly nonlinear waves. Fig. 2 shows two examples of 
two such solutions corresponding to the same input (marked by dots in the upper 
panel of Fig.l) along with their eigenvalue spectra. As expected, in both cases there 
is a continuum component filling densely the interval [—2, 2] on the real axis corre- 
sponding to propagation of linear waves 4 . In addition, isolated eigenvalues indicate 

4 Examining the eigenvalue spectra for different sizes M of the matrix, reveals that 
this continuum part of the spectra is mostly affected by truncation error. Typically 
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that the solution / undergoes an oscillatory instability while II is stable. The com- 
ponents of the corresponding eigenvectors are exponentially localized around the 
nonlinear portion of the chain. This is intuitively clear as only perturbations located 
there can destabilize the solutions. 

For comparison we also integrate the time-dependent DNLS equation setting as 
initial condition 0„(O) = u n and imposing the boundary conditions <j> n (t) = it„e~ lMt 
at the two edge sites n — —N p and n = N + N p to simulate the infinite system. As 
seen in Fig.3 the destabilization of solution / occurs by an exponential growth of 
<f> n on the central sites accompanied by an oscillatory behavior with the frequency 
prescribed by the stability spectrum (see the middle panel of Fig.3). Since this 
frequency is in the band of linear waves, this process is accompanied by emission of 
some radiation (traveling peaks in the upper panel of Fig.3) until the amplitude and 
frequency become large enough leading to a stable localized object (lower panel in 
Fig.3). This state is reminiscent of a nonlinear defect mode [18]. There is however 
an important difference as the localized mode is superimposed to a plane wave and 
that the overall evolution is quasi-periodic. 



4 Scattering of wave packets 

In this section we illustrate the consequences of the above results on the transmission 
of wave packets. In a nonlinear system where the superposition principle no longer 
holds, the connection between the two problems is not trivial. We solved numerically 
the time- dependent DNLS on a finite lattice \n\ < M with open boundary conditions, 
for the case of the dimer discussed in [21]. We take as initial condition a Gaussian 
wave packet (for no < 0) 



where fco > for no < (left-incoming packet) and where fco < for no > N (right- 
incoming packet). The upper panels of Fig. 4 display the evolution of two packets 
with the same I and opposite wave number fco impinging on the nonlinear dimer. 
The asymmetry of their propagation is manifest. In both cases, the packets are 
significantly distorted after scattering, and the emerging envelope may vary wildly. 
However, the Fourier analysis shows that they remain almost monochromatic at 
the incident wave number fco (lower panels of Fig. 4), with some small background 
amplitude radiation leaking throughout the lattice. 

To quantify the asymmetry of the scattering, and to compare with the above 
analysis, we measured the wave packet transmission coefficient as the ratio between 
the transmitted norm at the end tf in of the run divided by the initial one, namely 
(for n < 0) 



</>n(0) = -fexp 




+ ifcon 



(9) 




the numerical eigenvalues have a spurious imaginary part which is of order 1/M. 
This is not surprising since the corresponding eigenvectors are extended waves 
and thus more sensitive to boundary effects. 
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lattice site n 

Fig. 3. Upper panel: unstable evolution of the solution / depicted in Fig. 2; numer- 
ical integration of DNLS with initial condition corresponding to the stationary orbit 
with a small perturbation of size 10~ 5 to the site n = 1; Middle panel: evolution of 
Re<f>i: initial oscillations are at an angular frequency close to the imaginary part of 
the unstable eigenvalue, ReX = 1.316. At the later stage, a stable periodic oscillation 
sets in with frequency increases to 2.5 which is outside of the band of linear waves. 
Lower panel: snapshots of real (squares) and imaginary (diamonds) parts of central 
part the chain at t = 400. The upper curve is \4>„\ 2 showing the appearance of a 
localized excitation residing on the dimer. 
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Fig. 4. Numerical simulations of the propagation of Gaussian wave packets, Eq. (9) 
impinging on a DNLS dimer. Here Vo = -2.5, |feo| = 1.57, e = 0.05, M = 500, 
|J| 2 = 2.5, w — 50 and no = T250, respectively. Lower panels: Power spectra of the 
real part of <j> n at times t = (green dashed line) and t — 250 (red solid line). 



In Fig. 5 we show the transmission coefficient as a function of the packet intensity 
\I\ 2 for two widths of the initial packet. Comparison with Fig.l show that the region 
of wave amplitudes where the non-reciprocal behavior is maximal (that between the 
vertical dotted lines) corresponds qualitatively well to the range where packets with 
comparable intensities are transmitted more asymmetrically. 

Besides this, the wider packets show a highly irregular behavior in the sense that 
there is sensitive dependence on initial packets' parameters. This is somehow remi- 
niscent of chaotic scattering (see e.g. [31]). Empirically, we found that this behavior 
depends very much on the width uu: for small w the curves are much smoother. Qual- 
itatively speaking, this can be understood as follows. The dimer can be seen as an 
integrable, two-degrees of freedom system perturbed by a time-dependent force (the 
incoming wave packet) and subject to dissipation (the radiation towards the linear 
leads) [32]. If the packet is very narrow, the duration of the perturbation is short 
and one may argue that this just determines the initial condition for the dynamical 
system. The subsequent evolution will be almost regular and a slight change of the 
forcing will have a little effect. On the contrary, a wider packet will result in a more 
complex (possibly chaotic) dynamics yielding large changes in the transmission. 
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Fig. 5. Dependence of wave packet transmission coefficients on intensity for two 
different wave packet widths w. Other parameters as in the previous figure. 

5 Summary and conclusions 

We discussed the scattering problem for a linear Schrodinger chain with an embed- 
ded nonlinear, non-mirror symmetric dimer. The simplicity of the model allows to 
compute analytically a whole family of non reciprocal plane wave solutions, such as 
left- and right-incoming waves, with the same incident amplitude Ro and frequency 
[l, have sizeably different transmission coefficients t. 

We than addressed the question of dynamical stability for some specific values of 
Ro and [i. Solutions with a large enough t generically undergo oscillatory instabilities 
(see again Fig. 3). At the initial stage, the development of the instability occurs 
through growth and oscillation localized around the dimer sites, accompanied by 
emission of radiation towards infinity. Eventually, a stable extended state emerges 
coexisting with a nonlinear defect mode oscillating at a different frequency. This 
numerical finding suggests that the lattice may support exact stable quasiperiodic, 
nonreciprocal solutions. 

Numerical simulations demonstrated that equal Gaussian wave packets imping- 
ing on the dimer from the two opposite directions have indeed very different trans- 
mission coefficients t p . The packet amplitudes I for which such non-reciprocal be- 
havior is maximal qualitatively correspond to the range in which extended waves 
with the same wavenumber are transmitted more asymmetrically. Moreover, t p is 
sensitively dependent on the input parameters, like for instance the initial width w 
and amplitude I (see Fig. 5), a feature which is reminiscent of chaotic scattering. 



Left incoming 
Right incoming 
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